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(57) Abstract 



A non-invasive medical imaging tech- 
nique capable of evaluating, in situ, the oxy- 
genation state of body tissues (2) (e.g., by mea- 
suring the spectral properties of heme proteins) 
is described. The disclosed technique employs 
a multi-wavelength collimated source (1) and a 
collimated receiver (6) and performs a posi- 
tional and angular scan of the scattered radia- 
tion for each position of the incident beam. 
The resultant data is evaluated by employing 
imaging schemes which give differential 
weights to the contribution of various volume 
elements in the medium to the detector re- 
sponses at various frequencies. These measure- 
ments yield physiological information while 
being, for example, an indicator of physiologi- 
cal stress caused by disease or trauma. 
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A METHOD OF IMAGING A RANDOM MEDIUM 



10 



FIELD OF THE INVENTION 

This invention pertains to the detection and 
three- dimensional imaging of absorbing and/or 
scattering structures in complex random media by 
detecting scattered radiation emerging from the 
medium. 

BACKGROUND OF THE INVENTION 



The detection and identification of objects 
within a turbid medium has long challenged 

15 investigators from diverse technical disciplines. 
For example, at visible frequencies objects or 
vessels located in fog are rendered largely 
undetectable, thus inhibiting safe navigation. In 
material science methods have long been sought to 

20 remotely image faults or other inhomogenieties within 
optically transluscent or opaque materials. 

Objects buried in a random medium are 
rendered optically invisible because of the effects 
of scattering. Scattering causes light to experience 

25 multiple paths thereby obscuring phase information. 

In a general sense the detection or imaging 
of subsurface objects in a random medium requires the 
solution of an inverse problem. Specifically this 
requires an understanding of how the medium affects 

30 the propagation of light. Mathematically this can be 
described by a modulation transfer function. 

While prior art techniques have allowed 



35 
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" observers to detect the presence of an object within 

a turbid medium (see below) , none have allowed the 
observer to detect its depth or structure. In 
clinical medicine, for example, it would be useful to 
detect the presence, size, location, and depth of a 
5 suspected tumor. 

In particular, it would be especially 
desirable to perform such measurements in a manner 
which could also monitor the close association 
between oxidative metabolism and organ function. As 

10 described below, this relationship can be determined 

by measuring the oxygenation or redox-state of heme 
proteins in situ (e.g. hemoglobin, myoglobin, 
cytochrome oxidase) using optical transmission or 
reflectance techniques . 

15 An optical technique which could 

differentiate depth-of-f ield may be used to generate 
a 3-Dimensional image of body tissues which 
identifies the oxygenation state of those tissues. 
Such measurements would thus permit the 3-D imaging 

20 of the functional activity of the tissue. Indeed, 

the development of such a technique would be useful 
for other applications including detecting and 
imaging objects located within the atmosphere or, in 
oceanographic studies for the imaging of submerged 

25 vessels or other subsurface objects. 

Because of the strong scattering properties 
of tissues, it is often not practical to study their 
optical properties by a transmission measurement. On 
the other hand, by resorting to a reflectance or 

30 backscatter measurement the spectroscopic properties 
of tissue can be studied. 

Much interest in this technique has been 
generated in recent years because of its ability to 
monitor, in a continuous and non-invasive fashion, in 

35 situ, the oxygenation state of heme proteins. Such 
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° measurements are highly desirable because of the well 
established relationship between oxygen availability, 
oxidative metabolism and organ function. 

Using reflectance techniques, the effects 
that hypoxic hypoxemia, hypercapnia, carbon monoxide, 
5 and cyanide induced hypoxemia have on cerebral and 
cardiac energetics can be monitored and compared to 
simultaneous measurements of organ function. 

Such studies demonstrate the extreme 
sensitivity of organ energetics and function to 

10 minimal changes in tissue oxygen levels as measured 
by the oxygenation state of heme proteins. The 
latter findings are consistent with results obtained 
using phosphorus NMR. Recent clinical investigations 
employing NIR reflectance measurements have 

15 demonstrated that rapid changes in brain oxygenation, 
blood volume and energy state occurred in response to 
hyperoxia, moderate hypoxia and hypercapnia. These 
investigations emphasize the clinical usefulness of 
this technique. 

20 Although it has been recognized that optical 

studies of tissue (particularly in the reflectance 
mode) can potentially yield highly significant and 
useful information regarding organ function, such 
measurements are of little clinical value without the 

25 ability to characterize some degree of depth-of -field . 

In fact, depth-of-f ield differentiation is 
critical. Without a depth-of -field differentiation 
these techniques cannot even distinguish between a 
superficial bruise on the scalp and a more serious 

30 internal injury. 

The inability to differentiate 
depth-of-f ield in tissue by past optical measurement 
techniques has compelled the prior art to restrict 
measurements to the identification of the overall 

35 oxygenation state of the tissue (i.e., solely one 
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° dimensional information). For example, U.S. Patent 

No. 4,281,645 (Jobsis) discloses a spectrophotometry 
transillumination method for non-invasive monitoring 
of the metabolism of a body organ which performs this 
type of one-dimensional measurement. Jobsis also 
5 discussed . this prior art technique in Reflectance 

Spectrophotometrv and Non-Invasive, Infrared 
Monitoring of Cerebral and Myocardial 
Oxygen Sufficiency and Circulatory Parameters 
(1977). Such spectrophotometric techniques are not 

10 able to use the dispersion pattern generated by 

objects buried in the target medium to image the 
objects and establish their location within the 
surrounding target medium. 

Prior attempts have been made to locate and 

15 image objects located within a turbid medium by 

utilizing transillumination techniques, for example, 
diaphanography, where light is directed towards a 
target object. In transillumination the target 
object is detected by shadows in the light pattern 

20 emerging at the opposite side of the medium. 

Typical of such prior attempts include 
inventions disclosed by United Kingdom Patent 
Application Nos. 2,068,537A, 2,111,794, 2,154,731, 
and 2,092,856 and U.S. Patent No. 4,312,357. Such 

25 transillumination techniques only permit an observer 

to determine the two-dimensional outline of the 
target object. Thus its depth within the surrounding 
tissue and its three-dimensional structure remain 
unknown. Furthermore, these transillumination 

30 techniques require that a . specimen be relatively thin 
and that both sides of the specimen be accessible to 
the measurement device. 

U.S. Patent No. 4,555,179 (Langerholc et. 
al.) discloses a method and apparatus for the 

35 detection of objects in a scattering medium. The 
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° Langerholc technique employs a collimated light 

source which scans the medium. The reflected 
radiation is then analyzed to detect the presence of 
the target object. If the absorption characteristics 
of the 

5 target object were previously known, the observer 

could determine the object's depth within the 
medium. For clinical studies, such information can 
not be practically obtained. 

A further serious drawback of the Langerholc 

10 method is that, in practice, the method described can 
only locate objects relatively close to the surface 
of the medium being scanned. In addition, unlike the 
current invention, the measurement described involves 
integrating the backsca ttered signal over an area 

15 having a diameter equal to 2.5 times the thickness of 

the medium (i.e. the so-called "scattering zone"). 
This restriction renders such an approach useless for 
measurement of targets having limited geometries, 
i.e. the human body. Furthermore, by integrating 

20 over this area, specific information regarding 

variations in the position and angle dependent 
emerging flux is lost. 

On the other hand, the present invention 
recognizes the significance of performing a position 

25 and angle scan of the scattered light as essential in 

identifying regional variations in the absorptive and 
scattering properties of a turbid medium. 
Indeed, discussions of techniques such as 
Langerholc 1 s admit that they are only able, to detect 

30 objects along one dimension. These include wires or 
blood vessels. This is of limited value when the 
observer wishes to know the location, size, and shape 
of three-dimensional objects within a turbid medium. 

Bonner et al., in Model For Photon Migration 

35 in Tu rbid Biologic a l Medium (1987), recognized that 
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° analysis of positional information can be used to 
infer sub-surface properties of a random medium. 
However, Bonner does not disclose how such 
information may be used to generate an image of the 
target medium. 
5 While prior art techniques have allowed 

observers to detect the presence of an object within 
a turbid medium, none have allowed the observer to 
image its depth or structure. In many medical and 
other applications, such information is vital. 

10 The present invention describes an imaging 

technique which may' be used to study the brain and 
breast. The presence of tumors, cysts, hypoxic or 
infracted regions will be readily detectable. For 
obstetrical procedures, it will be possible to image 

15 the oxygenation state of an unborn fetus's brain, in 

utero , just prior to delivery. The progression of 
atherosclerosis on the delivery and utilization of 
oxygen by tissues (in particular, limbs) could be 
directly assessed. Because the technique of the 

20 invention is noninvasive and nondestructive and 

yields vital physiological information, it will also 
be useful in monitoring the response of body tissues 
to various therapies. In particular, the technique 
of the present invention may be used to monitor and 

25 evaluate the physiological status of burned patients, 
immuno-compromised patients or other patients 
restricted to isolation rooms. 

The inventive technique represents an 
accurate and reliable means to assess the functioning 

30 of transplanted organs as well as being a sensitive 
means to detect, in situ, impending organ rejection. 
The disclosed technique will be helpful in monitoring 
the physiological .status of excised donor organs for 
impending recipient implantation. Its use during 

35 anesthesia will permit, for the first time, direct in 
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u sit u monitoring of the delivery to and utilization of 
oxygen by the brain.' Such measurements in a hospital 
emergency room setting with unconscious patients 
would differentiate accurately and rapidly between 
patients suffering from carbon monoxide-, cyanide-, 
5 or other drug poisonings which interfere with oxygen 
utilization from patients having suffered a stroke or 
subtle cerebral trauma. 

Studies in marine environments using visible 
and/or near infrared (NIR) sources will permit the 
10 remote monitoring of the oxygenation status of plants 
and animals, i_n situ . Such measurements will 
indicate the physiological impact of environmental 
pollutants. Other marine studies may involve the 
imaging of turbulence at increasing ocean depths. 
15 The proposed method may also be used in search and 
rescue operations for the imaging of objects or 
victims buried in snow, ice or muddy waters. 

In a general sense, the above applications 
pertain to the imaging of objects which exist in the 
20 medium at the time of measurement. The present 

technique may also permit the determination of events 
which have previously occurred and thereby affected 
changes in the environment. This is particularly 
useful for the measurement of affects on marine life. 
25 Commercial applications of the proposed 

technique include the remote inspection of food 
products for spoilage or contamination by 
insecticides. As a remote imaging method, the 
proposed technique might aid aircraft or naval 
30 navigation in foggy atmospheres. Other commercial 
applications include the monitoring of various 
processes involved in the production of bulk 
industrial commodities at critical steps involving 
significant hazards to operators. The remote imaging 
35 of optically opaque objects which are not readily 
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° explored by a transmission measurement may also be 
imaged. Examples include the interrogation of low 
fault tolerant components of existing structures for 
the detection of fault lines or other fractures by 
the use of x-ray or particle beam sources. These may 
5 include nuclear reactor shielding, hulls of 

commercial or military aircraft, etc. 

Consequently, it is an object of this 
invention to provide a method for allowing an 
observer to accurately detect, three-dimensionally 
10 image and spectroscopically characterize target 
objects located within a turbid medium. 

Another object of the present invention is 
to employ radiation directed towards target objects 
located in a random medium and to detect radiation 
15 scattered from the medium to enable the observer to 

determine the object's depth, structure, absorptive 
and scattering properties within the turbid medium. 

A further object of this invention is to 
image reference objects in a non-invasive and 
20 non-destructive manner. 

Yet another object of the present invention 
is to provide a method whereby a physician may use 
the invention to aid in medical diagnosis. 

For example, a goal of this invention would 
25 be to measure the oxygenation state of body tissues 
and to display this information as a 

three-dimensional image to yield vital physiological 
information while still being a sensitive indicator 
of subtle physiological stress caused by disease or 
30 trauma. 

Additional objects and advantages of the 
invention will be set forth, in part, in the 
following description, will be obvious, in part, from 
this description, or may be learned from the practice 
35 of this invention. The objects and advantages of the 
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invention are realized and obtained by the processes 
and methods particularly pointed out in the following 
description and claims. 

SUMMARY OF THE INVENTION 

These and other objects are accomplished by 
the present invention which provides a method for 
detecting and imaging an object having different 
absorptive or scattering properties than that of its 
surrounding turbid medium. 

The inventive technique directs radiation 
into a medium towards, a suspected target or reference 
object located within the medium. The directed 
radiation may be multiple frequencies and may be 
generated by C.W. or ultra-fast pulse sources, (e.g. 
for time-of-f light studies). The invention 
recognizes that the analysis of position and angle 
dependence of scattered radiation inherently yields 
information which may be used for imaging; in 
particular, for the generation of three-dimensional 
images. 

In one aspect of the invention, by using 
reflectance spectrophotometry techniques, 
reflectance data from the reference object is 
analyzed. 

An aspect of the invention is the 
recognition that by directing a collimated beam of 
radiation into a turbid medium, radiation can be 
identified which .has propagated to successively 
increasing depths in the medium. This may be done by 
measuring positional dependent flux which emerges 
from the medium at increasing distances from the 
radiation source.. 

Positional dependent measurements may 
include measuring the angular dependence of emerging 
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° flux at increasing distance from the source. This 
off-axis measurement yields a "depth profile" of 
objects and other optical inhomogenieties within the 
medium. 

The invention recognizes that lateral 
5 detection, at greater distances from the incident 
source than previously contemplated, is a key to 
providing necessary imaging information. Whereas 
radiation back-scattered to a detector positioned 
laterally close to the incident beam is largely 

10 unable to detect deeply buried objects, such objects 
can be detected by receptors positioned at greater 
distance from the source. 

The invention uses data measured by a 
collimated receiver or receiver array which is 

15 directed towards the medium and located about the 
radiation source. The receiver array and light 
source are directed to multiple locations about the 
surface, thereby effecting a two-dimensional surface 
scan. 

20 Another aspect of the invention relates to 

the analysis of positional and angular dependent 
emerging flux. This data allows the recreation of a 
three-dimensional image of the medium. 

The invention recognizes that reflectance 

25 spectrophotometric data, processed correctly, 
inherently yields data which may be used for 
imaging. This is based on an observation that 
radiation undergoing multiple scattering at 
increasing depths within a turbid medium (such as 

30 living tissue) will emerge at increasing distance 
from the source. 

The invention recognizes that the proper 
analysis of data may involve a differential weighting 
of the targeted medium with each weighting function 

35 corresponding to a particular source detector 
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conf iguration . This differential weighting affords 
differential depth information regarding the 
scattering and absorptive properties of the medium. 

The invention allows the determination of 
the relative contribution of various volume elements 
5 within the medium to the detector response. 

Reconstruction of a three-dimensional image may be 
accomplished by the superposition of the relative 
contribution of each volume element, (i.e. voxel) to 
the detector response for each source-detector 
10 configuration. 

The invention may be utilized for various 
medical and non-medical applications. 

Therefore, in accordance with the present 
invention, a method for imaging a three dimensional 
15 object in a turbid medium is provided. The method 
comprises : 

A method for imaging a three-dimensional 
target object or other optical inhomogeniety in a 
turbid medium, comprising the steps of: 
20 providing a substantially collimated 

beam of radiation and directing the beam 
onto a surface of a turbid medium containing 
a target object; 

providing a substantially collimated 
25 receiver to receive the scattered radiation; 

performing positional scans of the 
scattered radiation over the surface of the 
turbid medium; 

applying the data obtained from each 
30 positional scan to form a three-dimensional 

image of the target object by: 



35 



a. determining attenuation of 
emerging radiation from said target 
medium relative to a model medium; 
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c b. determining a relative 

contribution of volume elements for 
each source-detector configuration; 
c* superimposing the relative 
contribution of volume elements for all 

5 source-detector configurations; 

d. repeating the above steps a to c 
for all desired frequencies to obtain a 
spectroscopic image . 



10 Another important aspect of the noviel 

technique is to perform a positional and angular scan 
of the radiation backscattered from the turbid medium 
wherein the target object is located. The technique 
employs a collimated light source and either a 

15 moveable collimated detector or a detector array, 

which permits the detection of a positional-dependent 
angular flux. The more differential the available 
information, the more which can be inferred about the 
properties of the backscattering medium. Moving the 

20 source relative to the detector yields additional 
information. 

BRIEF DESCRIPTION OF THE DRAWINGS 

25 The invention will be described in greater 

detail hereinafter by way of reference to the 
following drawings : 

Figure 1 is a graph illustrating the average 
maximum depth of penetration of emerging light as a 

30 function of radial distance from the source. 

Figure 2 is a graph illustrating the depth 
distribution of emerging photons at different 
distances from the source as a function of maximum 
depth of penetration of emerging. 

35 Figure 3 is a graph illustrating the 
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reduction in light intensity at increasing distances 
from the source in a multilayer agar medium. 

Figure 4 is a graph illustrating a radial 
difference plot of the effect of a buried absorber on 
the intensity of emerging light from a Monte Carlo 
calculation . 

Figure 5 is a graph illustrating a radial 
difference plot of the effect of a buried absorber on 
the intensity of emerging light from experimental 
results. 

Figure 6 is an abbreviated illustration of 
the interaction of the technique of the present 
invention with a turbid medium. 

Figure 7A is an abbreviated illustrative 
rendition of the contour of a relative weight 
function level utilized by the method of the 
preferred embodiment of the invention. 

Figure 7B is an abbreviated rendition of the 
superposition of the contour of figure 7A. 

Figure 7C is an abbreviated rendition of the 
image reconstruction performed by the technique of 
the preferred embodiment of the invention. 

Figure 8 is an abbreviated flow-chart of the 
post-receiver information processing of the method of 
the preferred embodiment of the present invention. 

DETAI LED DESCRIPTION OF THE PREFERRED EMBODIMENT 

The measurement and analysis routine of the 
present invention may be readily appli.ed to the 
imaging of any scattering medium. The inventive 
technique involves inherently remote imaging which 
may be directed to not only living tissue, but 
oceans, atmospheres or other turbid media. 

The radiation employed by the invention is 
not restricted to visual radiation, but may be 
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radiation of much higher or mucli lower wavelength 
including the entire spectral range and other types 
of radiation, e.g., x-ray, audible radiation or 
wavelengths used in radar. 

An aspect of the invention introduces a 
5 method for imaging a three-dimensional target object 
or other optical inhomogeniety in a turbid medium, 
comprising the steps of: 

providing a substantially collimated 
beam of radiation and directing the beam 
10 onto a surface of a turbid medium containing 

a target object; 

providing a substantially collimated 
receiver to receive the scattered radiation; 
performing positional scans of the 
15 scattered radiation over the entire surface 

of the turbid medium; 

applying the data obtained from each 
positional scan to form a three-dimensional 
image of the target object by: 

20 

a. determining attenuation of emerging 
radiation from said target medium 
relative to a model medium; 

b. determining a relative contribution 
25 of volume elements for each 

source-detector configuration; 

c. superimposing the relative 
contribution of volume elements for all 
source-detector configurations; 

30 d. repeating the above steps a to c 

for all desired frequencies to obtain a 
spectroscopic image. 

The inventor has recently characterized the 
35 three dimensional dispersion of light in a random 
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° medium of isotropic scatterers with a buried absorber 
using a Monte Carlo simulation routine. These 
results, described below, serve to illustrate the 
importance of considering positional and angular 
dependent emerging flux as a means to selectively 
5 interrogate increasing depths in a random medium. 

The mean value + standard deviation 
(expressed in terms of mean free paths (mfp) of the 
average maximum depth to which the emerging light 
penetrates) are shown in Figure 1 as a function of 

10 radial distance from the source. Here, the 

simulation involved launching one million photons 
from a point collimated source and allowing each 
photon to undergo up to 20,000 collisions. The mean 
free path was equal to unity. 

15 The dispersion of photons in such a medium 

as described by Figure 1 is simulated by generating 
three random numbers from a pseudo-random number 
generator. These numbers represent the angles & 
(theta), and^(phi), and the distance R. O is the 

20 angle between the vector represented by a scattered 
photon and the Z axis. ^ is the angle between the x 
axis and the vector sum of the x and y components of 
the vector represented by the scattered photon. R is 
the free path length of the photon to the next 

25 collision within the dispersive medium. 

Each time a new scattering vector is to be 
computed for a photon within the medium three new 
random numbers are needed; one each for 0, ^ and R. 
The random number generator produces . a number between 

30 0 and 1. The uniform scattering model assumes that 

photons uniformly scattered from the center of a unit 
sphere must produce a uniform density of photons that 
reemerge from the surface of the sphere. Phi must 
therefore be uniformly distributed between 0 and Iff 

35 radians (360 degrees). A random number, n, from the 



WO 89/12223 



PCT/US89/02511 



-16- 

random number generator will have a value uniformly 
distributed between 0 and 1. Multiplied by 2*7/^ This 
will result in a random number with an equal 
probability of falling anywhere within the range of 0 
to 2<7f This is used for the value of jz( 
5 Theta must be distributed between 0 and if 

such that its cosine has an equal probability of 
falling anywhere within the range of -1 to +1. 
Accordingly, 

# 

10 cos &= 1 - 2 * n (i) 

The free path R is distributed according to: 



L. 



P(L) = 1 e 



where is the mean free path for the medium. 

20 

Given random number n, set 

L = L * Log -r 1 — (3) 
m ^ 1-n v 1 

25 Equation (3) will yield a random variable L 

distributed as specified by the function P(L) in 
equation (2) . 

There is a near linear relationship between 
the average maximum depth of penetration of the 

30 emerging light and the distance from the source. 
Thus, light emerging farther from the source has 
preferentially penetrated to a greater depth than 
light emerging near the source. 

The fraction of incident photons which have 

35 propagated to their maximum depth prior to emerging 
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° at a distance of 1 and 10 mfp from the source as a 

function of maximum depth of penetration is shown in 
Figure 2. This fraction was determined by standard 
Monte Carlo techniques. At a distance of 1 mfp from 
the source, approximately 60% of the emerging light 
5 has penetrated to a maximum depth of between 0-1 mfp 
into the medium followed by an exponential decline at 
greater maximum depths. In contrast, at a distance 
of 10 mfp from the source, most of the emerging light 
has propagated to a maximum depth of 4-8 mfp into the 
10 medium. The conditions of the simulation were the 
same as described in Figure 1. 

Experimental results shown in Figure 3 show 
the relative intensity of backscattered light, at 
increasing distances from the source of multilayer 
15 turbid agar media with biological absorbers, compared 
to an agar medium with no absorbers. 

The medium with no absorbers was prepared by 
adding microscopic uniform latex beads (0.202 urn 
diameter) to liquified clear agar and allowed to gel, 
20 forming a 6mm thick layer. Agar medium containing 

latex beads and intact, heat treated red blood cells 
was similarly prepared. The intensity of 
backscattered light was then measured for the agar 
medium containing three layers of no absorbers and 
25 compared to agar media having a middle layer of 

intact red blood cells and a middle and bottom layer 
of intact, heat treated red blood cells respectively. 

The turbidity of the top agar layer rendered 
the lower absorbing layers undetectable to the eye 
30 when viewed in room light. However, as shown in 

Figure 3, these layers become increasingly detectable 
when measured by a collimated source (He-Ne laser) 
and a collimated receiver positioned laterally at 
increasing distances from the source. 
35 Further evidence substantiating the merits 
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of an off-axis measurement are shown in Figure 4. 
Therein, results show a Monte Carlo calculation 
demonstrating the effect of a buried absorber on the 
intensity of emerging light as displayed by a radial 
difference plot. This simulation involved launching 
of 1 million photons from a point collimated source 
at coordinates 50,50 with each photon undergoing a 
maximum of 400 collisions. Photons emerging through 
an acceptance angle of 120-150° with respect to the 
direction of launching were detected at each point 
along the surface of 100x100 mfp grid. The absorber 
formed a 2x2x2 mfp cube and was positioned 8 mfp from 
the source and 4 mfp below the surface. The plot 
shown in Figure 4 was smoothed by calculating the 
moving average over a 3x3 mfp area. It is. apparent 
that the impact of the buried absorber is best seen 
by evaluating intensity data at increasing distances 
from the source. 

A corresponding evaluation of experimental 
results demonstrating the effect of a buried absorber 
on the intensity of emerging light as displayed by a 
radial difference plot is shown in Figure 5. A two 
dimensional surface scan was performed using a 
collimated receiver tilted 30° from normal. 

In Figure 5, Z-coordinate values correspond 
to the intensity difference, in percent, between one 
side of the source and its radial opposite. Values 
for x and y axis are indicated in thousandths of an 
inch from the source. The medium consisted of 
uniform latex beads (0.202 urn diameter) in water in 
which a black plastic bead (0.1" diameter) was 
suspended at a depth (0.75") sufficient to render it 
undetectable from the surface when viewed in room 
light. 

The lateral displacement of the absorber 
from the source is indicated by the open circle. The 
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solid circle indicates the x, y coordinates of the 
calculated "center-of-mass" of the intensity 
differences on one-half of the plain. Analysis of 
the change in position of the calculated 
"center-of-mass" for several beam positions correctly 
5 determined the position of the absorber to within 10% 
of its actual location. 

The invention provides a unique technique 
for imaging complex media by initially assigning to 
every point (locations and direction) in the medium a 
10 weight given by the product of the flux (direct 

solution) and the expected contribution of a photon 
at the point to the detector response (adjoint or 
backward solution) for the medium without the 
absorber. 

15 If the absorption is small, the effect of an 

absorber at that point is proportional to the 
weight. The effect of all the absorbers is a linear 
combination of such contributions. When the 
absorption is not small, the effect is no longer 

20 linear, but the linear combination will serve as a 

first approximation. To obtain the weight functions, 
one must compute both the direct and adjoint fluxes, 
or alternatively, directly compute the contribution 
flux. 

25 Determination of these weight functions for 

specified media are relatively straightforward, 
although for unknown heterogeneous media such as body 
tissues the solution is more complex (see below). 
However, the anatomy of the body is well defined. 

30 Large veins near the surface are readily apparent and 
the location of other major arteries and veins deep 
within tissues are also well known. In the context 
of an imaging scheme, these major components of the 
vascular tree can viewed as optical landmarks which 

35 would aid in the calibration of the weighing 
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° functions. 

Clearly, the more accurately the 
contribution of various regions of tissues can be 
assigned to the detected 

signal, the greater the image resolution and 
5 spectroscopic accuracy. The detection of these 

optical landmarks together with the measurement of 
angular dependence of emerging flux (as a function of 
distance from source and position of source relative 
to the media) will provide differential information 

10 concerning the optical properties of tissues. The 

more such differential information is available, the 
more detailed the inferrences about the properties of 
the medium which can be made. A time-to-flight 
analysis could also be performed in conjunction with 

15 the weight function. 

In a simple consideration of the weight 
functions, one might decide that when the weights are 
normalized to run from 0 to 100, an absorber at any 
point with a weight of less than 5 cannot be 

20 detected. One keeps only the volume within this 

contour and assumes that any detectable absorber must 
lie inside. A given source-detector reading is 
observed relative to the reading expected without an 
absorber. Any absorber must lie in the corresponding 

25 volume. The absorber, if any, must then lie in the 
volume complementary to the union of all of these 
excluded volumes. 

This incremental exclusion process provides 
an image of the absorber. For absorbers near the 

30 surface the image is sharp, because of the 

contribution by single scattering. However, the 
resolution will degrade for succeeding depths. One 
may repeat this procedure continually, excluding 
volumes outside the contour surface of weight 10, 20 

35 or 50. Of course, each tightening of the acceptable 
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° volume tightens any image obtained at the possible 
expense of missing it altogether . 

In its preferred embodiment, the device of 
the present invention would detect light scattered 
from an optical inhomogeneity , such as a tumor, 
5 located within an optically dispersive medium such as 
normal human tissue. The incident radiation in this 
embodiment would be near infrared (NIR) light with a 
wavelength of 750 - 1000 nm. 

The instrumentation of the present 
10 invention, in a preferred embodiment, is controlled 
by a microcomputer and have, as major components, an 
NIR tunable dye laser, a precision linear translator 
having 5 degrees of freedom, two NIR-sensit ive 
photomultiplier tubes, phase-locked amplifiers and 
15 electronic components to detect the light signal. 

Functionally, the instrument operates 
similar to a double-beam scanning spectrophotometer. 
A beam splitter will direct a portion of the incident 
light from the dye laser to an NIR-sensit ive PMT 
20 (reference signal), while the remaining light will be 
directed to the sample. Scattered light will then be 
detected by a second NIR-sensi tive PMT housed within 
a collimator (sample signal) . 

The latter device may be attached to an 
25 overhead precision linear translator having 5 degrees 
of freedom of movement (x, y, z, phi, theta). 
Signals from the two PMT's will then be directed to 
the amplifiers and the signal ratio determined. The 
intensity of the laser, its output wavelength, 
30 positioning of the linear translators, and 

acquisition and storage of the measured light signal 
may be controlled by a microcomputer. 

While the above discussion has emphasized 
the ability to image variations in the absorptive 
35 properties of random media, the technique is equally 
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w capable of imaging variations in the intensity of 
backscattered light caused by differences in the 
scattering cross section of the medium unrelated to 
its absorptive properties. Differentiation between 
these two effects can be made by examining the 
5 spectroscopic properties of the medium, i.e., 
comparing results at different wavelengths. 

In a second preferred embodiment of the 
invention, a measurement scheme employs a source of 
collimated radiation in the near-infrared region 

10 (NIR) of the spectrum. 

The source directs radiation towards a 
target medium (i.e., body tissue). Radiation 
propagating through the medium will experience 
multiple scattering prior to emerging. 

15 A collimated receiver or receiver array 

detects the emerging radiation. The receiver array 
may be translated to perform a position and angle 
scan of the radiation. These scans are performed 
about the collimated source which is also translated. 

20 The present invention characterizes the 

relative contribution of emerging flux from each 
volume element in the medium. The characterization 
considers the depth of the volume element and 
distance and angle of the receiver from the source. 

25 As described below, the technique relates the 

information obtained from actual measurements of 
positional dependent flux to determine the relative 
contribution of each volume element to the detector 
response. This information is used to reconstruct a 

30 3-dimensional image. 

In the shown preferred embodiments, the 
geometry of the source-detector configuration is 
arranged to perform a reflectance (i.e., backscatter) 
measurement. However, additional embodiments are 

35 envisioned in which the geometry of the 
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° source/detector configuration allows transmission and 
other measurement . schemes to be performed. 

In another preferred embodiment the source 
may utilize either a continuous or ultra-fast pulse 
laser which may operate at various frequencies. The 
5 use of ultra-fast sources with a streak-camera as a 
receiver allows time-of-f light measurements of 
intensity. The analysis of time-of-f light data is . 
similar to that for a continuous source. However, 
the advantage of the time-of-f light measurement would 

10 be to reduce the volume of medium contributing most 
to the detector response, thereby yielding a higher 
resolution image . 

Figure 6 depicts the interaction of the 
technique of another preferred embodiment with a 

15 random turbid medium. A collimated incident beam of 
radiation is directed towards a turbid medium 2. In 
this preferred embodiment incident beam 1 may be a 
multif requency source in the near infrared (NIR) 
region of the spectrum. This is preferred for 

20 spectroscopic studies of body tissues because at 

these frequencies light is known to penetrate deeply 
into the tissue medium. Additionally, a 
mult if requency measurement can be utilized to 
identify the relative oxygenation state of heme 

25 proteins. 

Radiation propagating through the medium 3 
experiences multiple collisions 4 prior to emerging 
from the medium 2. 

Emerging radiation 5 will be detected by a 
30 collimated receiver array 6 positioned about the 

source 1 thereby effecting a surface scan. The angle 
of the receiver 6 relative to the source 1 may also 
be varied 7. 

Both the source 1 and receiver array 6 may 
35 also be translated about the surface 8 in order to 
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° perform an overall surface scan and thereby greatly 
improve subsequent image resolution. 

The source 1 and the receiver array 6 are 
collimated, and caused to undergo their respective 
surface and angular scans to enhance subsequent image 
5 resolution. 

In this preferred embodiment of the 
invention, receiver arrays 6 are positioned about the 
source 1 to reduce scanning time. However, other 
embodiments are envisioned which engage either a 

10 single collimated receiver as described above or a 

single receiver array which performs multiple scans. 
Additionally, in a further embodiment of the 
invention multiple sources may be utilized to mimick 
the translation of a single source. 

15 In this preferred embodiment, the receiver 

array 6 may be caused to contact the medium's surface 
or may be position at a distance from the target 
medium 2 when detecting emerging radiation 5. In one 
embodiment of the invention the receiver arrays are 

20 cooled CCD elements, however the receiver may be any 
type of radiation intensity detection mechanism. The 
detected radiation intensity is transformed into an 
electrical signal and transmitted 9 to a data 
processing unit 10. The data processing unit 10 

25 transforms the radiation intensity signal into a 
3-dimensional image of the target medium. . 

Figures 7 and 8 illustrate the processing of 
information by this preferred embodiment of the 
invention. The flow-chart of Figure 8 indicates how 

30 the detected intensity information by the receiver 
array (indicatory numeral 6 in Figure 6) is applied 
as an electrical signal to the data processing scheme. 

In the data processing scheme 3 the measured 
intensity of emerging radiation. is compared to the 

35 expected intensity thereby determining an attenuation 
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° coefficient of a corresponding scattering volume. 

The expected intensity is determined from calculation 
or actual measurements of a defined model medium as a 
function of position and angle of the receiver 
relative to the source. These calculations can be 
5 readily performed using a Monte Carlo simulation. 

The intensity attenuation coefficient is 
then applied to a weight function to determine the 
relative contribution of volume elements of the 
target medium to the detector response for each 
10 source-detector configuration. 

Image reconstruction for an unknown medium 
is based on weight functions calculated for a defined 
medium. A given location in the space of the target 
medium has a preferential contribution or weight to a 
15 detector response for radiation propagating in the 
random medium. The contribution of all points in 
space to the detector response is determined by 
calculating the product of flux through a given voxel 
and the probability the photons passing through this 
20 voxel will contribute to the detector response. 

By this model, the effect of weak absorption 
in the medium on the detector response is assumed to 
be a linear combination of the product of the 
absorption cross-section and corresponding weight 
25 function for the affected voxels. A linear 

approximation is considered for simplicity of 
explanation of the preferred embodiment. However, in 
practice a non-linear approximation will yield a more 
accurate spectroscopic analysis and an image of 
30 greater resolution. 

A non-linear model considers an interative 
scheme in which the absorption in a given voxel 
causes a reduction in flux through the neighboring 
volume elements. This scheme applies corrections to 
35 the weight functions calculated in model media to 



WO 89/12223 



PCI7US89/02511 



-26- 

° more accurately estimate the true contribution of all 
points in space to the detector response for the 
target medium. 

Calculation of the weight functions can be 
readily performed using the Monte-Carlo method. The 
5 result yields a three-dimensional map of the relative 
contributions of all points in space to the detector 
response for a given source-detector configuration. 

in practice, to the extent the values of the 
weight functions for a defined medium differ 

10 significantly from the actual values for the unknown 
medium, the resultant image will appear distorted. 
Provided the distortion is not too severe, 
recognizable optical landmarks (i.e. major vessels) 
can be located within the tissue. The location of 

15 these landmarks can be compared with the known 
anatomy (as determined, for example, by an MRI 
image) , thereby serving as a basis for image 
enhancement and deconvolut ion. 

The resolution of the images obtained will 

20 be greatest near the surface and decline in concert 
with the uncertainty associated with the values of 
the weight functions at greater depths. A 
2-dimensional depiction of a relative weight function 
contour map of the target medium is illustrated by 

25 Figure 7A. Therein can be seen a source 201 directed 
to the surface of a turbid target medium 202. The 
receiver 203 detects emerging radiation scattered by 
the medium 202. Inside the medium 202, radiation 
propropagating through each one of the volume 

30 elements (voxels) 204 will all contribute to the 
detector 203 response. The envelope 205 is a 
schematic illustration of a representative contour. 
This contour is an arbitrarily designated fraction of 
weight functions of those voxels contributing most to 

35 the detector response. 
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As depicted by the illustrations shown in 
Figure 7, an image of the subsurface structure in a 
random medium can be obtained by calculating the 
product of the attenuation coefficient (as determined 
by comparing the measured intensity to an expected 
5 intensity from a model medium) and the corresponding 
weight function for each voxel for each 
source-detector configuration. The sum of these 
products corresponds to the superposition of all 
voxels for all source-detector configurations. This 

10 calculation, as a function of all source positions, 
yields a three dimensional map of the relative 
absorption cross-sections over all space (Figure 7B) . 

The image may then be reconstructed to view 
a 3-dimensional depiction of the spectroscopic 

15 properties of the target medium. The spectroscopic 
image may be generated by comparing each of the 
reconstructed images obtained at each of the various 
frequencies. This is done using standard 
spectroscopic techniques. 

20 A pictoral representation of a spectroscopic 

image may be seen in figure 7C. This reconstructed 
image may be displayed graphically on paper or on a 
video screen or any pictorial representation known in 
the art . . 

25 While the preferred embodiments and 

techniques of the invention have been presented in 
detail, modifications and adaptations of such 
techniques and embodiments will be apparent to those 
skilled in the art. However, it is to be expressly 

30 understood that such modifications and adaptations 
are within the spirit and scope of the present 
invention as set forth in the following claims. 



35 
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CLAIMS 

I Claim: 

1. A method for imaging a 
three-dimensional target object or other optical 
5 inhomogeneity in a turbid medium, comprising the 
steps of: 

providing a substantially collimated 
beam of radiation and directing the beam 
onto a surface of a turbid medium containing 
10 a target object; 

providing a substantially collimated 
receiver to receive the scattered radiation; 

performing positional scans of the 
scattered radiation over the surface of the 
15 turbid medium; 

applying the data obtained from each 
positional scan to form a three-dimensional 
image of the target object by: 

a. determining attenuation of 

20 emerging radiation from said target 

medium relative to a model medium; 

b. determining a relative 
contribution of volume elements for 
each source-detector configuration; 

25 c. superimposing the relative 

contribution of volume elements for all 
source-detector configurations; 
d. repeating the above steps a to c 
for all desired frequencies to obtain a 

30 spectroscopic image; 

displaying said spectroscopic image. 



2. The method of Claim 1, wherein said 
35 beam of radiation is near infrared radiation. 
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3. The method of Claim 1, wherein said 
collimated receiver is a receiver array. 

4. The method of Claim 3, wherein said 
receiver array includes CCD detectors. 

5. The method of Claim 1, wherein said 
receiver is positioned for a reflectance measurement, 



6. The method of Claim 1, wherein said 

10 receiver is positioned for a transmission measurement 

7. The method of Claim 1, wherein said 
positional scan includes having said collimated beam 
of radiation be laterally scanned, 

15 

8. .The method of Claim 1, wherein said 
positional scan includes a lateral translation of 
said collimated receiver. 

20 9 - The method of Claim 1, wherein said 

positional scan includes an angular translation of 
said collimated receiver to provide differing angular 
information about the detector response. 

25 10. The method of Claim 1, wherein the data 

is applied to a data processing apparatus, 

11. The method of Claim 1, wiierein said 
collimated receiver detects information regarding 

30 intensity of said scattered radiation. 

12. The method of Claim 1,. wherein 
calculation of attenuation is performed by comparing 
measured values to expected values from. a model 

35 medium calculated using Monte Carlo techniques. 
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° 13. The method of Claim 12, wherein said 

weight function includes a calculation of a product 
of flux through a given voxel and an expected 
contribution of said voxel to a detector response. 

5 14. The method of Claim 1, wherein said 

attenuation step includes determining the ratio of 
actual data obtained from each positional scan to a 
calculated expected detection from said model medium. 



10 15. The method of Claim 1, wherein said 

determining relative contribution step includes 
multiplying said determined attenuation by a value of 
the weight function of each volume element for each 
source-detector conf iguaration under consideration. 



15 



20 



16. The method of Claim 1, wherein said 
superimposing step determines a three dimensional 
image of said target at a single frequency 
measurement . 

17. The method of Claim 1, wherein said 
repeating step includes superimposing of all results 
of said superimposing steps. 



25 18. The method of Claim 1, wherein said 

model medium is simulated by assigning a weight to 
every point in the medium, wherein said weight is a 
product of the flux and expected contribution of a 
photon at a point, to the detector response for a 

30 medium without an absorber. 

19. The method of Claim 1, wherein said 
imaging is used to evaluate, in situ, body tissue 
oxygenation state. 
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20. The method of Claim 1, wherein said 
imaging is employed to measure spectroscopic 
properties of heme proteins, in situ. 
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